VORTICITY AND ENTROPY PRODUCTION IN TILTED SZEKERES 
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We analyze the properties of the tilted Szekeres spacetime, i.e. the version of such spacetime as 
seen by a congruence of observers with respect to which the fluid is moving. The imperfect fluid 
and the kinematical variables associated to the four-velocity of the fluid assigned by tilted observers 
are studied in detail. As it happens for the case of the Lemaitre-Tolman-Bondi spacetime, the 
fluid evolves nonreversibly (with nonvanishing entropy production) and is nongeodesic. However 
unlike the latter case, the tilted observer detects vorticity in the congruence of the fluid world lines. 
Also, as for the nontilted congruence the magnetic part of the Weyl tensor vanishes, reinforcing the 
nonradiative character of this kind of spacetime. Possible physical implications of these results are 
discussed. 
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I. INTRODUCTION 

In any physical theory, part of the description of a spe- 
cific scenario is tightly related to the congruence of ob- 
servers carrying out the study [l| . In the case of general 
relativity, the existing arbritariness in the choice of the 
four velocity in terms of which the energy-momentum 
tensor is split leads to a variety of physical interpreta- 
tions for the source of a given spacetime (see and 
references therein). 

As it has been shown in the references above, when the 
two possible interpretations of a given spacetime corre- 
spond to two congruences of observers related by a local 
Lorentz boost, both the general properties of the source 
and the kinematical properties of the congruence would 
be different. 

Particularly enlightening is the case of tilted Lemaitre- 
Tolman-Bondi (LTB) spacetime whose fluid source is 
nongeodesic and evolves nonreversibly, unlike the non- 
tilted version (see [l?} for details). 

An important (though obvious) remark is in order at 
this point: we are assuming general relativity to be valid. 
This implies that special relativity is always valid at the 
local level. Therefore we can always perform a Lorentz 
boost (locally), thereby ensuring the existence of tilted 
observers. 

In this work we endeavour to study in detail tilted 
Szekeres spacetimes [HI, [l|| . The motivations to under- 
take such a task are many. 



Szekeres dust models have no Killing vectors and 
therefore represent an interesting generalization of LTB 
spacetimes [2(| |2l|. They have been extensively used as 
cosmological models (see [22l . [23j j and references therein) 
and in the discussion of the gravitational arrow of time 

Physical and topological properties of Szekeres space- 
times, as well as the problem of structure formation and 
other cosmological related issues have been considered in 

MM. 
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A renewed interest in Szekeres spacetimes has ap- 
peared, in relation with the potential role of density in- 
homogeneity to explain the observed cosmic acceleration, 
without invoking dark energy (see [35r!42j and references 
therein) . 

As we shall see below, a heat flux term will appear in 
the energy-momentum tensor of the tilted congruence. 
As it happens for LTB pj} , such a heat flux vector would 
be necessarily associated to a "truly" (i.e. entropy pro- 
ducing) dissipative phenomenon, even in the absence of 
bulk and shear viscosity. 

We shall see that as in the nontilted version, the mag- 
netic part of the Weyl tensor vanishes, implying that 
the super-Poynting vector asociated to the Bel-Robinson 
tensor vanishes too, which suggests the absence of grav- 
itational radiation. Thus confirming the nonradiative 
character of Szekeres spacetime already stressed in (20j . 

Besides the dissipative (irreversible) evolution of the 
tilted version, and the nongeodesic character of the fluid, 
another important difference with respect to the non- 
tilted version appears, namely: the vorticity of the fluid 
lines defined by the tilted four-velocity is now nonvanish- 
ing. As it can be expected, this vorticity has important 
consequences on the physical interpretation of the model. 
We shall comment on this issue in the last section. 
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II. SZEKERES SPACETIME 

Before considering the tilted version of Szekeres space- 
time we find it useful to provide a very brief review of 
this spacetime in its standard, nontilted, version. 

Since t he p ublication of the first two original papers by 
Szekeres [l8l . [l9} this metric has been presented in many 
different forms. Here we shall closely follow (with slight 
changes in notation) the formulation given in [26]. 

There are two classes of Szekeres spacetimes, which in 
Szekeres notation correspond to ft ^ and /3' = 0, here 
we shall deal only with the first class, denoted as LT-type 
Szekeres metric in [26J. 

Thus, the line element is given by: 



ds 2 = -dt 2 



where a prime denotes derivative with respect to r, R = 
R(t,r), e = ±1,0 and / = f(r) > — e is an arbitrary 
function of r. We number the coordinates x° = t, x 1 = 
r, x 2 = p, x 3 = q. 

The function E is given by 



E(r,P, q) 



S 
2 



P 



P 



S 



q-Q 

S 



(2) 



where S = S(r), P — P(?') and Q = Q(r) are arbitrary 
functions. 

From Einstein equations it follows that R satisfies the 
equation 



R 2 



2M 



+ /, 



(3) 



where a dot denotes derivative with respect to t, and 
M = M(r) is an arbitrary function. From the above 
equation it follows that 



R 



-M 



R 2 



(4) 



from where the meaning of M as an effective gravitational 
mass becomes evident. 

As in the LTB case there are three subfamilies of so- 
lutions depending on the value of /, namely; hyperbolic 
(/ > 0), parabolic (/ = 0) and elliptic (/ < 0). 

On the other hand, the parameter e determines how the 
2-surfaces, r = constant, foliate the 3D spatial sections 
of constant t. These may be spheres (e = +1), planes 
(e = 0) or quasispheres (e = —1). 

For the "spherical" case e = +1, we have the follow- 
ing relationship between coordinates q,p and "spherical 
angle" coordintes (0,4>). 



(P-P) 



cot - cos < 
2 



(Q-Q) J ■ 

— = cot - sin < 

S 2 



(5) 



(6) 



The LTB spacetime is recovered in the case E' = 0, 
e = +l. 

For the nontilted (comoving) congruence the source of 
Szekeres spacetime is pure dust, therefore the energy- 
momentum tensor is: 

T ll u = ^{t,r,p,q)v ll v v , (7) 

where vector i> M (V'u^ = —1) is just 

^ = (1,0,0,0), (8) 

and p, denotes the energy momentum density as measured 
by comoving observers. Then, from the (t, t) component 
of the Einstein equations we have 

2(M' - MIE'/E) 



8lTfl R 2 (R' — RE' jE) ' 

The two nonvanishing kinematical variables are: 
The expansion scalar 



= ^ 



R' E — RE' R 



R'E-RE 1 W 
and the shear tensor, which can be written as: 

1 



a n u n» 



-h„ 



where h^ v = + v^v v , a^a^ = |cr 2 and 



(n^rift — 1), with 



g <^>^ ,3 
' R'E - RE' ' ' 



E RR' - RR' 



R \ R'E-RE' J ' 
Next, the electric part of the Weyl tensor 



can be written as 



where 



R \ii> £ ( n^n v r^iiv 



R 

e = 3- + 4tt/2. 

R 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



The magnetic parts of the Weyl and the Riemann ten- 
sors vanish identically. 

The fact that both the shear tensor and the electric 
part of the Weyl tensor have only one nonvanishing inde- 
pendent component, as well as the vanishing of the mag- 
netic part of the Weyl tensor illustrates how close this 
spacetime is to the spherical symmetry. This similarity 
is further reinforced by the fact that the specific forms 
([TTj) and ([15)) are typical of spherical symmetry. Thus 
the qualification of "quasispherical" assigned by Szekeres 
himself to his solution [lj| looks well justified, in spite of 
the fact that it has no Killing vectors. 

We shall next consider the tilted version of Szekeres 
spacetime. 
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III. TILTED SZEKERES SPACETIME 

In order to obtain a tilted congruence, let us perform 
a Lorentz boost from the Locally comoving Minkowskian 
frame (associated to w M ) to the Locally Minkowskian 
frame with respect to which a fluid element has "radial" 
(in the r direction) velocity lj. For simplicity we shall 
consider a boost only in the "radial" direction and lj to 
be a function of t and r alone. 

The corresponding tilted congruence is characterized 
by the four-velocity vector field 



1 



#(e + /) 1/2 

;i - lj 2 ) 1 / 2 ' R'E - RE' (1 - to 2 ) 1 / 2 



,0,0 



(17) 

For the tilted observer, the matter distribution de- 
scribed by ([7]) will be that given by 



where use has been made of 
1 



(1 - a; 2 ) 1 / 



1/2 Va 



(1-lj 2 ) 1 / 



T75 JV a . 



(25) 



Alternatively, it may be useful to define two auxiliary 
variables (P r and Pj_) as: 



P r = N a N T a t 



P ± = K a K^T at 



(26) 



where K a is a unit spacelike vector (orthogonal to V a ) 
in the p direction. 
Then we may write 



n Q/3 = n(N a Np - -hap), 



(27) 



Tfiu = ftVpVv + PV + n ^ + Qn v <> + V »9" ( 18 ) 

where fx is the energy density, P the isotropic pressure, 
g M the heat flux vector and II Ml , the anisotropic pressure 
tensor, that this observer measures (of course the projec- 
tor tensor is now defined in terms of V 1 instead of 

The above is the canonical, algebraic decomposition 
of a second order symmetric tensor with respect to unit 
timclikc vector, which has the above indicated standard 
physical meaning when T a p is the energy-momentum ten- 
sor describing some energy distribution, and the four- 
velocity assigned by certain observer. 

It is immediate to see that: 



~^V a -T aP V^ 



(19) 



with 



n = P r - Pj 



p = 



Pr + 2P ± 



(28) 



from where the physical meaning of P r and Pj_ becomes 
evident. 

Using the above definitions, the relationships linking 
tilted and nontilted variables, besides ([21]) and (|22|) . read 



Pi 



0: 



3P = IT = P r 



/MjJ 



(29) 



We have next calculated the kinematical variables for 
the tilted congruence. These were obtained with MAPLE 
13. 

The four-acceleration is now nonvanishing, and its 
components are: 



p = -h af) T af! , n„^^(T F -pg. (20) 

A quick calculation taking into account the above, ren- 
ders: 



j 2 EE! - E'R 



(1 - lj 2 ) 2 R'E - RE 1 1-lj 2 R'E -RE' {1-lj 2 ) 2 '' 

(30) 



1-LJ 2 



P = 



)ILJ 
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(21) 



u R'E- RE' ' 



(31) 



q a =qN a , q = - 



fjjjj 

1 - LJ 2 



(22) 



2 

a = 



w 2 A E 



1 -lj 2 R 



a 3 = 



lj 2 t e 



1 -lj 2 R 



(32) 



N a 



" ^ + /) 1/2 I o 0^ 

(1-w 2 ) 1 / 2 ' R'E -RE' {1-lj 2 ) 1 / 2 ' ' J' 

(23) 



U aP = -^^(N a Np - \h a p), (24) 

1 — LJ O 



where subscripts p, q denote derivative with respect to 
these coordinates and 



_ E' p E - E'E p 



E'E — E' E n 

r = -4r^ — (33) 



R'E - RE' ' R'E - RE 

The expansion scalar now becomes 
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1 ERR' + R(2R'E - 3RE') 



U10J 
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u' E^/TTJ 
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I 



(34) 



Next, for the vorticity tensor 
1, , 1 



a. 



2 {V KV - Vu.,,) + - ( av y v - a u V^) 



we find the following nonvanishing components 

w 2 A R u 2 T R 



fin 



12 



n 13 - — 



(l_ w a)3/2 E ' 
uj R 

~( 1 _ w 2)3/2 E 

uj R 



fin 



(1-uj 2 ) 3 / 2 E 



(35) 



,(36) 



y/e + f{E'E-tfE p ), (37) 



(1-w 2 ) 3 / 2 E 
whereas for the vorticity vector 



-^7Tf(E'E-E'E q ), (38) 



L 7' 

we find 



a, = n^upa Vc.fi V°g va g pl3 ; 



r /0123 



fi = Oi = o, 



fi 2 = 



1 - uj 2 E 



fi a 



-wA i? 



1 -cj 2 £ 



(39) 

(40) 
(41) 



Finally, for the shear tensor the following nonvanishing 
components appear, 



°00 = 77 



3 1 -u 
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(42) 
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(43) 
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(44) 
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522(700, (45) 



(46) 



where 



(7/ 



(1-0,2)3/2 



i? 



£ RR' - R'R 



LJUJ Uj' E\Jt + f 

(1 - w 2 ) 1 / 2 R R'E - RE' ' (1 - a; 2 ) 3 / 2 + (1 - uj 2 fl 2 R'E - RE' ' 



+ 



(47) 



It is a simple matter to check that the shear tensor may 
be written as: 



Q /3 = (Xi(N a N - -h a p) + a H {N a K p + NpK a ) 



+ (Tni{N a Lp + N p L a ), 



where 



luT 



on 



2(1 -w 2 )' 



2(1 -uj 2 )' 



(48) 
(49) 



and L a is a unit spacelike vectors in the q direction. 

Observe that comparing (|24|) with (j48|) it follows that 
the anisotropy described by II Q ^ cannot be associated 
exclusively to shear viscosity. 

It is worth stressing that the appearance of vorticity 
is unavoidable for a generic tilted Szekeres spacetime. 
Indeed, vanishing of vorticity implies because of (|4"Tj) 



A 



0. 



(50) 



r 



producing 



Eg 



(51) 



where ip is an arbitrary function of its arguments. Then 
using © in (f5Tj) we obtain 



p - P{ r ) = q) - Q{ r )^iv, q) 



(52) 



which implies P(r) — Q(r) = 0, where the regularity con- 
dition P(0) = Q(0) = has been used. Conditions (|5D|) 
can also be satisfied for some axially symmetric Szkcrcs 
models [43| . 

Thus for a generic Szekeres spacetime there is a non- 
vanishing vorticity according to the tilted observer. Fur- 
thermore, as it follows from (|49l) the shear tensor has 
the form (|48[) . implying three independent components 
instead of one for the nontilted observer. 
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It is worth mentioning that the magnetic part of the now has nonvanishing components which are associated 



Weyl tensor vanishes, like for the nontilted case. However 
the magnetic part of the Riemann tensor 



to the the dissipative flux q a . They are: 



Z„ 



1 



a0 R va afi V»V° 



(53) 



^23 



Z: 



R 



1 



32 



1 - w 2 E 2 R'E - RE 



- [E(RR' + 2RR') + E'(2R 2 + RR)}. 



(54) 



Finally, the electric part of the Weyl tensor 
now has the following components: 



2 J 1 
Eoo — - o £ > 



2,1-lo 1 



En 



31-w 2 



"g 322 ^, 



E; 



33 



E, 



where 



£ 



E R R — RR 1 
R R'E - RE' ' 



(55) 

/<7ii£,(56) 
(57) 

(58) 



From the above it follows that we may write the electric 
part of the Weyl tensor as: 



E a0 =E(N a N - ^h a0 ), 



(59) 



I 

which is similar to the expression (fl~5|) for the nontilted 
case. 



IV. TRANSPORT EQUATION AND ENTROPY 
PRODUCTION 



In order to evaluate the possibility of entropy produc- 
tion we need a transport equation for the heat conduc- 
tion. 

Due to the well known pathologies of Eckart [44| and 
Landau 45[ approaches we shall resort to causal dissi- 
pative theories. Reasons for doing that have been exten- 
sively discussed in recent years (see \U> 5!) and references 
therein) . 

Thus for a second order phenomcnological theory for 
dissipative fluids we obtain from Gibbs equation and con- 
servation equations (see [H, [H[ for details) : 



TS C , 



(60) 



where S a is the entropy four-current, T is temperature 
and Pi = -Sf, where r and k denote thermal relaxation 
time and thermal conductivity coefficient respectively. 
Then from the second law of thermodynamics 



S" > 0, 



(61) 



the following transport equation is obtained (see [5a , l58j 
for details): 

Th^Vq^ + q a = -nh afi (Tp + Tap) 

For simplicity, in the above equations we have not 
included thermodynamic viscous/heat coupling coeffi- 
cients, neither have we included couplings of heat flux to 
the vorticity. Also, because of the comment below (|4"9")l 



and for simplicity we have omitted the contributions of 
shear viscosity ( see 53[ for details). 

Let us first consider the situation within the context 
of the standard (Eckart-Landau) irreversible thermody- 
namics, in which case r = 0. Then after simple manipu- 
lations we obtain 



5 a 
■n 



1. 



(63) 



where C x denotes the Lie derivative with respect to 



the vector field x" 



V" 
T 



and T. 



a/3 



dis. 



From the above it is evident that if \ defines a conformal 
Killing vector (CKV), i.e. 



for an arbitrary function $, then 
S" = 0, 



(64) 



(65) 
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a well known result [oOj . 

Since a generic Szekeres spacetime does not admit 
CKV [sij, then the heat flux q a has to be associated 
to an irreversible process. 

In the case of the causal thermodynamics the following 
expression emerges from (|60[) 

S- a = -\Tf s .H x 9a, ~ \ {^f) ■ (66) 

Due to the presence of different phenomenological pa- 
rameters in different terms in (|66[) , it becomes evident 
that for a generic Szekeres spacetime S a a =^ 0. 

We recall that in the above expression we have ne- 
glected terms involving couplings of heat flux to the vor- 
ticity. These couplings give rise to terms of the form 
K,T"fxQ a pqP , where 71 is the thermodynamic coupling co- 
efficient, which have to be added to the right hand side 
of ([62]) ■ Also, the omitted shear viscosity contribution in 
([50)1 would give rise to a term of the form 2r]^2^ t ^ a ul P> ^ 
which has to be added to the transport equation for the 
shear viscosity, where 72 is a coupling coefficient, r\ de- 
notes the shear viscosity coefficient and <> is the spatial 
tracefree part of the tensor. 

Thus, as in the LTB case [l7|, tilted observers in Szek- 
eres spacetime detect a real (entropy producing) dissi- 
pative process, while for the nontilted observer the evo- 
lution proceeds adiabatically. In this case however it is 
worth noticing that vorticity contributes to entropy pro- 
duction too. 

Also, as in the LTB case we might speculate that the 
origin of such important difference in the pictures be- 
tween both sets of observers may be found in the fact that 
while the fluid is geodesic for nontilted observers, it is not 
for tilted ones. Then invoking the equivalence between 
collisional terms and force terms in the Boltzmann equa- 
tion established in [62j , it could be possible that the force 
term associated with the nonvanishing four-acceleration 
plays the role of a collisional term, leading to entropy 
production. 

V. CONCLUSIONS 

We have analyzed the hydrodynamic and thermody- 
namic properties of the tilted Szekeres spacetime and 
compared them with the nontilted case. 

Two main differences emerge from our study: Firstly, 
the fluid as seen by the tilted observer evolves irreversibly 
in contrast with the picture described by the nontilted 
observers according to which the system evolves adia- 
batically, furthermore in the former case the fluid is no 
longer geodesic. These differences between both set of 
observers also appear in the FRW and LTB spacetimes. 

Secondly and more interesting, the fluid lines as seen 
by the tilted observers have nonvanishing vorticity. As 
mentioned before this is true for any generic Szekeres 
spacetime. 



From a cosmological point of view this might have in- 
teresting consecuences. Indeed, the possibility of rotat- 
ing universes (besides the well known case of the Godel 
model) and eventual restrictions on the magnitude of 
the associated vorticity have been studied in the past 
(see [6314661 ] and references therein) . Also different mech- 
anisms for vorticity generation in a cosmological back- 
ground have been presented (see [67147H and references 
therein). Here we have seen that the appearance of vor- 
ticity may be also an observer related phenomenon. 

We would like to conclude with the following remarks: 

• We have considered the simplest possible case of 
tilted congruence, it is obvious that a boost in all 
possible directions would produce a much richer 
physical picture from the point of view of tilted ob- 
server. In the same line of arguments we have con- 
sidered the simplest dissipative scenario and have 
not considered the possible presence of an electro- 
magnetic field. The reasons for doing that are on 
the one hand the cumbersome resulting expressions 
and on the other, the fact that even at this level of 
simplicity very interesting physical results emerge. 

• It is worth noticing that, as it follows from (|29l) . 
if jl has a good physical behaviour, so do have the 
physical tilted variables (at least for small values of 
W ). 

• In spite of the differences, the tilted version con- 
serves some of the "structural" properties of the 
nontilted one, e.g. the magnetic part of the Weyl 
tenssor vanishes and the electric part has only one 
nonvanishing independent component. 

• For a generic Szekeres spacetime with no symme- 
tries and non-barotropic equation of state, it is not 
obvious that the Gibbbs framework would be appli- 
cable (see 72] ) . This point deserves more attention. 
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